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We present a quark model calculation of the charmonium spectrum with self energy corrections 
due to the coupling to the meson-meson continuum. The bare masses used in the calculation are 
computed within the relativized quark model by Godfrey and Isgur. The strong decay widths of 
3S, 2P, ID and 2D cc states are also calculated, in order to set the values of the 3 Po pair-creation 
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model's parameters we use to compute the vertex functions of the loop integrals. Finally, the 
nature of the X(3872) resonance is analyzed and the main possibilities (cc state or DD* molecule) 
are discussed. According to our results, the A(3872) is compatible with the meson Xci(2P), with 
, J PC = and is thus interpreted as a cc core plus higher Fock components due to the coupling 

D ■ to the meson-meson continuum. These J PC = 1 ++ quantum numbers are in agreement with the 

experimental results found by the LHCb collaboration. In our view, the A (3872) 's mass is lower 
, than the quark model's predictions because of self energy shifts. 
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I. INTRODUCTION 



a,- 

The quark model (QM), in all its possible reformulations |l!4ll|. can properly describe several properties of the 
hadrons, such as the spectrum and the magnetic moments, but it neglects continuum coupling effects. Indeed since 
the earliest days of hadron spectroscopy, it has been recognized that properties of levels can be strongly influenced 
[ — ■ by nearby channels The presence of these higher Fock components in meson and baryon wave functions are 

lf~) ' predicted by the QCD and must have an effect on the QM similar to that of unquenching lattice QCD calculations. 
, In particular, these continuum coupling effects can contribute, through a self-energy term, to a shift in the hadron 
' masses, as already shown by several authors in the baryon (l3l - [l9| and meson [20l - [3l| sectors. 



CNJ , Interest in loop corrections in the meson sector 22j was triggered after the discovery of the narrow charmed-strange 
mesons D* (2317) + and D s i(2460) + [Hj], since their surprisingly low masses could be explained by this type of 
effect. In the 80 's, Tornqvist et al. [2fj | studied heavy cc and bb quarkonium within the unitarized quark model and 
calculated the mass shifts and mixing induced by DD, D*D* , ... loop diagrams, using the 3 Po decay model [H| for 
hadron vertex functions. 

Barnes and Swanson [27| computed the mass shifts of charmonium IS, 2S and IP resonances due to DD, DD*, 
D*D* , D S D S , D S D* S , D* S D* S loops. The authors evaluated the coupling between the valence component and the 



continuum component by using the 3 Pq model 3^|, with Gaussian meson wave functions. 

Danilkin and Simonov analyzed the mass shifts of charmonium A^Si (TV = 1, 2, 3) [28[ and 2 3,1 Pj [29[ states, using 
the mechanism of channel coupling via decay products. The authors applied the Weinberg eigenvalue method [1^] to 
multichannel problems, considering DD, DD* and D*D* decay channels. 

Eichten et al. [2l| evaluated the influence of open-charm channels on charmonium properties, such as strong decay 
widths and self energies. The authors revisited the properties of charmonium levels, using the Cornell coupled-channel 
model [36| to assess departures from the single channel potential-model expectations. 

Hwang and Kim [22| calculated the mass shift of D*j(2317) due to coupled channel effects, within the Cornell 
coupled-channel model of Rcf. j36|. According to them, the measured mass of this meson, being 160 MeV lower 
than the corresponding estimation of Ref. [H , appears surprisingly low and can only be explained by coupled channel 
effects. 

The loop corrections can be relevant to the study of the AT (3872) meson [37j], whose nature has not yet been 
understood. Indeed, there are currently two possible interpretations for the meson: a weakly-bound 1 ++ DD* 
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molecule [Ifl [29|, [3(1 [H, HH or a cc state liol . l4l| . with 1 ++ or 2 h quantum numbers. For a summary of theoretical 
interpretations of the X(3872), see Ref. [421 ] . 

In the last few years, interest in heavy mesonphysics has increased enormously, as has the number of collaborations 
devoted to the topic. In particular, BaBar [H,|4J], Belle [45j], CDF [46| and DO have already provided many interesting 
results; moreover, all four detectors at LHC (Alice, Atlas, CMS and LHCb) have the capacity to study charmonia 
and bottomonia and have already produced some results, such as the discovery of a Xb(3P) system [47j |. There are 
also approved proposals for new experiments, such as Belle II (48[. 

The calculation presented in this article is the first attempt to calculate in a systematic way the spectrum of 
charmonia within a quark model, including loop corrections, and makes it possible to perform a comparison with the 
already existing and the future experimental data. Something similar has already been done for bottomonia in Refs. 

Our results for the spectrum of charmonia are fitted to the experimental data [49( , so that the calculated masses of 
the mesons of interest are the sum of a bare energy term, computed within the relativized QM by Godfrey and Isgur 
[1] , and a self energy correction, computed thanks to the formalism of the unquenched quark model (UCQM) [HI, [H(| . 
In our UCQM calculation, we consider as intermediate states a complete set of accessible SUf(4) ground-state (i.e. 
IS) mesons. IS intermediate states, being at lower energies than P-wave and -D-wave intermediate meson states, give 
the main contribution to the self energies of the charmonium states that we are going to study. 

Furthermore, we present some results for the strong decay widths of charmonium 35, 2P, ID and 2D states, 
calculated within a modified version of the 3 Po pair-creation model [13|. This is done in order to set the values of the 
3 Po model's parameters we use to compute the vertex function of the UCQM [see Eq. (J25J]. 

Finally, we use our results for the cc spectrum to discuss the nature of the X(3872) resonance. Specifically, we 
analyze the interpretation of this meson as a cc state with 1 ++ or 2 h quantum numbers. According to our results, 
the X(3872) is compatible with the meson Xd(2 3 Pi), with J PC = 1++. 



II. FORMALISM 
A. Self energies 

The Hamiltonian we consider, 

H = H + V , (1) 

is the sum of an "unperturbed" part, Hq, acting only in the bare meson space, and of a second part, V, which can 
couple a meson state to a continuum made up of meson- meson intermediate states. 

The dispersive equation, resulting from a nonrclativistic Schrodingcr equation, can be written as 

W = EfA^f, (2a) 

BC JO — J^bc 

where the bare energy E a satisfies: 

M a = E a + £(P Q ) ■ (2b) 

M a in Eq. ()2b[) is the physical mass of the meson A, with self energy T,(E a ). In Eq. (|2a[) one has to take the 
contributions from various channels PC's into account. A channel PC is a meson-meson intermediate state, with 
relative momentum q and quantum numbers Jb c and £ coupled to the total angular momentum of the meson A. The 
matrix element V a .bc of Eq. (|2a|) results from the coupling, due to the operator V, between the intermediate state 
PC and the unperturbed quark-antiquark wave function of the meson A; Ebc = Eb + E c is the total energy of the 
channel PC, calculated in the rest frame of A. Finally, if the bare energy of the meson A, i.e. E a , is greater than the 
threshold Eb c , the self energy of Eq. (|2"aj) contains poles, and is a complex number [see Eq. (|15p]. 

Since the physics of the dynamics depends on the matrix elements V r aj b c (g), one has to choose a precise form for the 
transition operator, V, which is responsible for the creation of qq pairs: our choice is that of the unquenched quark 
model of Ref. [3l|, so a 3 P model. 



B. Unquenched quark model 

In the unquenched quark model [3ll . |50T ] the effects of qq sea pairs are introduced explicitly into the quark model 
(QM) through a QCD-inspircd 3 Pq pair-creation mechanism. This approach, which is a generalization of the unitarized 
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quark model by Tornqvist and Zcnczykowski jl3| . was motivated by later work by Isgur and coworkers on the flux- 
tube breaking model. They showed that the QM emerges as the adiabatic limit of the flux-tube model to which the 
effects of qq pair creation can be added as a perturbation (52|. Therefore, our approach is based on a QM to which 
the quark- antiquark pairs with vacuum quantum numbers are added perturbatively. The pair-creation mechanism 
is inserted at the quark level and the one-loop diagrams are computed by summing over the possible intermediate 
states. 

Under these assumptions, the meson wave function is made up of a zeroth order quark-antiquark configuration plus 
a sum over all the possible higher Fock components due to the creation of 3 Pq quark-antiquark pairs. To leading order 
in pair creation, the meson wave function is given by 



4>a) = M 



A) + d v\ BCqiJ) 

BCIJ J 

(BCgtJ I Tt I A) ' 
E a — Eb — E c 



(3) 



where represents the 3 Pq quark-antiquark pair creation operator [53}, A is the meson, B and C are the intermediate 
virtual mesons, and E a , Ef, = yj M? + q 2 and E c = */ M 2 + q 2 are their respective energies, q and I the relative radial 
momentum and orbital angular momentum of B and C , and J is the total angular momentum, with J = J& + J c + 1. 



The 3 Pq quark-antiquark pair-creation operator of Eq. ([3]) can be written as [52 



T f = -370 Jdp 3 dpid(p 3 + pi) C34 -F34 e r; 
[X34 x yi (p 3 -p 4 )}^ bl{p 3 ) 4(p 4 ) , 



-P4) 2 /6 



(4) 



where b 3 (p3) and d\(J>4) are the creation operators for a quark and an antiquark with momenta p 3 and P4, respectively. 
The quark and antiquark pair is characterized by a color singlet wave function C34, a flavor singlet wave function 
-F34, a spin triplet wave function X34 with spin S = 1 and a solid spherical harmonic IVi(p3 — Pi) that indicates that 
the quark and antiquark are in a relative P wave. Since the operator creates a pair of constituent quarks with an 
actual size, the pair creation point has to be smeared out by a gaussian factor, whose width r q was determined from 
meson decays to be in the range 0.25 — 0.35 fm [T7L . 

The pair-creation strength 70 is a dimcnsionless constant, fitted to the strong decay widths of cc states (see Sec. 
IIII Al for details). The matrix elements of the pair-creation operator were derived in explicit form in the harmonic 
oscillator basis as in Ref. [53| . using standard Jacobi coordinates. The meson wave functions have good flavor 
symmetry and depend on a single oscillator parameter a, which, according to the literature [13, EH, [5(|, is taken to 
be a = 0.50 GeV. 





B 



FIG. 1: Two diagrams can contribute to the process A — >• BC. a and qi stand for the various initial (i = 1 — 4) and final 
(i = 5 — 8) quarks or antiquarks, respectively. Picture from Ref. [31]]. APS copyright. 



In the UCQM, the coupling V a ,bc between the continuum channel BC and the unperturbed wave function of the 
meson A can be calculated as 



V aM {q) =J2{BCqU\T^\A) 



(5) 
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In general, two different diagrams can contribute to the transition matrix element {BCq£J\ \A) (see Fig. [TJ: in the 
first one, the quark in A ends up in B, while in the second one it ends up in C. In the majority of cases, one of these 
two diagrams vanishes; however, for some matrix elements, both must be taken into account [3l| . as for example, 
this is the case of the coupling r\ c — > J / ^ J / ^ , where the initial |cc) state is coupled to the final state \cc;cc) and the 
created pair is a cc one. 

Finally, the expression for the self energy of the meson A, Eq. (|2aj). can be re- written as 



BCiJ 



C. Godfrey and Isgur's relativized quark model 



There is a hug e number of studies on meson spectroscopy, based on different pictures for mesons; they include qq 
mesons 0, 0, meson- meson molecules [H, 0, US, [39l I5p - l60| . tetraquarks [HI, 0, EH and quarkonium hybrids 
[63Tj65j and references can be found in reviewpapers like (66l |. 

The relativized QM by Godfrey and Isgur |2J is a potential model for qq meson spectroscopy. This model assumes 
a relativistic dispersion relation for the quark kinetic energy, a QCD- motivated running coupling constant a s (r), a 
flavor-dependent potential smearing parameter er, and replaces factors of quark mass with quark kinetic energy. 

The Hamiltonian of the model [2J is given by 



H = 



^conf + ^hyp 



so 



(7) 



where mi and m,2 are the masses of the constituent quark and antiquark inside the meson, q is their relative momentum 
(with conjugate coordinate r), V^ on f, ^Jryp an< ^ are ^ ne confining, hyperfine and spin-orbit potentials, respectively. 
The confining potential, 



V, 



conf 



3 3 L 

- c + - or 

4 4 



a s {r) 



(8) 



contains a constant, c, a linear confining term and a Coulomb-like interaction, depending on the renormalized running 
coupling constant of QCD, a s (r) (for more details see Ref. Q); moreover, one has: 



(qq\ Fx • F 2 \qq) 



(9) 



The hyperfine interaction is written as 



hyp 



The spin-orbit potential 0], 



I 3 Si r S 2 -r 



f S 1 -S 3 6 3 (r) 



(10) 



Vso - Vs ,cm + V S0)tp , 



(11) 



is the sum of two contributions, where 



V s , 



r- 5 V rrti m j J 

(A + A.' 

\ mi rrij 



L Ft ■ Fj 



(12a) 



is the color-magnetic term and 



_ldH™f (St sA - 



(12b) 



is the Thomas-prcccssion term. 



5 



III. RESULTS 



A. Strong decay widths 



In this section, we show our results for the strong decay widths of 35, 2P, ID and 2D charmonium states above 
the DD threshold (see Table Hill). 

The decay widths are calculated within the 3 Po model [I?], EH, [5(| as 

T A ^ B c =$A^Bc(qo)^2\(BCq lJ\Tl \A)\ 2 . (13) 

$a->bc(?o) is the standard relativistic phase space factor [55l[56j. 

^ E b (q )E c (q ) 

®A^>bc = 2nq — , (14) 

Ma 

depending on the relative momentum qo between B and C and on the energies of the two intermediate state mesons, 
E b = y/Mg + ql and E c = y/M? + ql (for the values of M b and M c , see Table HJ). The operator T T inside the 3 P 



State Mass [GeV] Source 



D 1.867 [49] 

D*(2007) 2.009 [49J 

D s 1.969 [49J 

D* s 2.112 [49J 



TABLE I: Masses of open charm mesons used in the calculations. 



amplitudes (BCqo IJ\T^ \ A) is that of Eq. Q, which also contains the quark form factor of Refs. [52|, [54|. The 
introduction of this quark form factor, which is just a Gaussian function in the relative momentum between the quark 
and the antiquark of the created pair, in the 3 Po model transition operator determines slightly different values for 
the model parameters (see Table HI)) . Specifically, the value of the pair-creation strength 70, which is fitted to the 
reproduction of the experimental strong decay widths of Table IIII1 is greater than that which would be obtained in 
the standard 3 Pq model [EflUfl], i.e. 70 = 0.4. 



Parameter Value 



70 0.510 

a 0.500 GeV 
r q 0.335 fm 

m„ 0.330 GeV 

m s 0.550 GeV 

m c 1.50 GeV 



TABLE II: Parameters of 3 P the model. 

Another difference between our calculation and those of Refs. [H, [5(| is the substitution of the pair-creation 
strength 70 with the effective strength 7^ of App. [B] The introduction of this effective mechanism suppresses those 
diagrams in which a heavy qq pair is created. More details on this mechanism can be found in Refs. [23L l3lj. 

Finally, the results of our calculation, obtained with the values of the model parameters of Table [TT1 are reported 
in Table Mil This set of parameters is also used in the self energy calculation of Sec. MI Bl in order to compute the 
vertices (BCq£J\ \A) of Eqs. © and |(T5J. 
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State DD DD* D* D* D S D S D S D* D*D* Total Exp. 



this 1 So) 


- 


38.8 


52.3 


- 


- 


91.1 


- 


* (4040) (3 3 Si) 


0.2 


37.2 


39.6 


3.3 


- 


80.3 


80 ± 10 


h c (2 1 Pi) 




64.6 








64.6 




Xco(2 3 P ) 


97.7 










97.7 




Xc2(2 3 P 2 ) 


27.2 


9.8 








37.0 




*(3770)(l 3 Di) 


27.7 










27.7 


27.2 ± 1.0 


cc(l 3 D 3 ) 
cc(2 1 J D 2 ) 


1.7 










1.7 






62.7 


46.4 




8.8 


117.9 




#(4160) (2 3 Di) 


11.2 


0.4 


39.4 


2.1 


5.6 


58.7 


103 ±8 


cc(2 3 D 2 ) 




43.5 


49.3 




11.3 


104.1 




cc(2 3 D 3 ) 


17.2 


58.3 


48.1 


3.6 


2.6 


129.8 





TABLE III: Strong decay widths (in MeV) for 3S, 2P, ID and 2D charmonium states. The values of the model parameters 
are given in Table HU The symbol - in the table means that a certain decay is forbidden by selection rules or that the decay 
cannot take place because it is below threshold. 



B. Bare energy calculation within the relativized quark model. Self energies of cc states 

The relativized QM Q, which is described in Sec. Ill C[ is here used to compute the bare energies of the cc states 
that we need in the self energy calculation. In our study, we computed the bare energies E a 's of Eq. (|2b|) as the 



m c = 1.562 GeV b = 0.1477 GeV 2 a" = 0.600 

A = 0.200 GeV c = 0.069 GeV a = 1.463 GeV 

s = 2.437 e c = -0.2500 e t = 0.0300 

e so( v) = -0.0314 e so(s) = 0.0637 



TABLE IV: Values of Godfrey and Isgur's model parameters, obtained by fitting the results of Eq. I|2b[) to the experimental 
data [11]. 

eigenvalues of Eq. (|7J|. At variance with QM calculations, such as that of Ref. 0, we did not fit th e eig envalues of Eq. 
d?]) to the experimental data |4jj. In our case, the quantities fitted to the spectrum of charmonia [491] are the masses 
M a 's of Eq. (|2b[) and therefore the fitting procedure is an iterative one. Our resulting values for the parameters of 
Godfrey and Isgur's model are shown in Table IIVI 

Once the values of the bare energies are known, it is possible to calculate the self energies H(E a )'s of IS, 2S, IP, 2P 
and ID cc states through Eq. If the bare energy of the meson A is above the threshold BC, i.e. E a > Mf, + M c , 
the contribution to the self energy due to the meson-meson channel BC is computed as 

Z(Ea]BC) 

- TQ + m c A \(BCqiJ\T^ \A)\ 2 (15) 

\ {BCqU\T^\A)\ 2 } ^ 



2-ni 



Eh c —E a 



where the symbol V represents the principal part integral, which can be computed numerically, and 

2ni j gE ^ Ec I (BCq I J\ I A) 1 2 1 is the imaginary part of the self energy, related to the decay width by: 

L * > E bc =E a 

T A ^ BC = Im[E(E a ;BC)} . (16) 

Finally, the results of our UCQM calculation, obtained with the set of parameters of Tables HT1 and IIVI and with the 
effective pair-creation strength 7^ of App. [51 arc shown in Table fVl and Fig. [2] 
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M 

(GeV) 

4.0' 

3.5- 

3.0 

2.5' 

2.0' 

1.5' 



X(3872)_ 



o-+ t— 1+- fr+ 1++ 2"+ 2— 3~ J 



PC 



FIG. 2: Comparison between the calculated masses (black lines) of 15, 2S, IP, 2P and ID charmonium states via Eq. (|2b[l 
and the experimental ones [4!| (blue boxes). The experimental mass of the X(3872) is indicated by a yellow box. The new 
values of the parameters of Godfrey and Isgur's model are taken from Table [TV] 



State 




DD 


DD* 


D*D* 


D S D S D S D* S D* s D* a 


VcVc 


Vc J/<f J/*J/* E(K) 


E a 


M a 










DD* 






D s D* a 


















^(l'So) 


o-+ 




-34 


-31 




-8 


-8 






-2 


-83 


3062 


2979 


2980 


J/*(l 3 Sl) 

Vc^ 1 So) 




-8 


-27 


-41 


-2 


-6 


-10 




-2 




-96 


3233 


3137 


3097 


0"+ 




-52 


-41 




-9 


-8 






-1 


-111 


3699 


3588 


3637 


#(2 3 Si) 
hc(l x Pi) 




-18 


-42 


-54 


-2 


-7 


-10 




-1 




-134 


3774 


3640 


3686 


1+- 




-59 


-48 




-11 


-10 




-2 




-130 


3631 


3501 


3525 


Xco(l 3 Po) 


0++ 


-31 




-72 


-4 




-15 







-3 


-125 


3555 


3430 


3415 


Xcl(l 3 A) 


1++ 




-54 


-53 




-9 


-11 






-2 


-129 


3623 


3494 


3511 


Xc2(l 3 P 2 ) 

h c {2 1 Pi) 


2++ 


-17 


-40 


-57 


-3 


-8 


-10 







-2 


-137 


3664 


3527 


3556 


1+- 




-55 


-76 




-12 


-8 




-1 




-152 


4029 


3877 




Xco(2 3 Po) 


0++ 


-23 




-86 


-1 




-13 







-1 


-124 


3987 


3863 




Xci(2 3 A) 


1++ 




-30 


-66 




-11 


-9 






-1 


-117 


4025 


3908 


3872 


Xc2(2 3 P 2 ) 

cc{l 1 D 2 ) 


2++ 


-2 


-42 


-54 


-4 


-8 


-10 







-1 


-121 


4053 


3932 


3927 


2" + 




-99 


-62 




-12 


-10 






-1 


-184 


3925 


3741 




*(3770)(l 3 L>i 


) 1- 


-11 


-40 


-84 


-4 


-2 


-16 









-157 


3907 


3750 


3775 


cc(l 3 D 2 ) 


2— 




-106 


-61 




-11 


-11 




-1 




-190 


3926 


3736 




cc(l 3 D 3 ) 


3— 


-25 


-49 


-88 


-4 


-8 


-10 




-1 




-185 


3936 


3751 





TABLE V: Self energies, E(_B a ) (in MeV, see column 12), for charmonium states due to coupling to the meson-meson continuum, 
calculated with the effective pair-creation strength of Eq. (|B1[) and the values of the UCQM parameters of Table [Til Columns 
3-11 show the contributions to ~E(E a ) from various channels BC, such as DD, DD* and so on. In column 13 are reported the 
values of the bare energies, E a , calculated within the relativized QM 0], with the values of the model parameters of Table HVl 
In column 14 are reported the theoretical estimations M a of the masses of the cc states, which are the sum of the self energies 
E(_E a ) and the bare energies E a . Finally, in column 15 are reported the experimental values of the masses of the cc states, as 
from the PDG H. 



C. Nature of the X(3872) resonance 



The quark structure of the X(3872) resonance, observed for the first time by the Belle Collaboration in the decay 
of the B meson 37 1 and then confirmed by CDF [67j], DO [68| and BABAR (69|, still remains an open puzzle. Indeed, 
at the moment, there are two possible interpretations for the meson: a weakly bound 1 ++ molecule [25l. l29l. [30L l38l. |39| 
or a charmonium state, with 1 ++ or 2 ^ quantum numbers 14 ill . For sure, we can only say that the decay channels 
where it was observed imply 1 ++ or 2 ^ quantum numbers [701 ]. while the other hypotheses are excluded by more 
than 3cr [7l|. It is thus necessary, in order to study properties of the X(3872) such as the decay modes, to make an 
assumption regarding its quark structure that is compatible with the quantum numbers 1 ++ or 2 
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The first and easiest possibility is to consider the A(3872) as a cc state [4l|. In this case, the X(3872) would 
correspond to a 2 3 P 1 resonance [x c i(2-P), J PC = 1 ++ ] or to a \ l D 2 (J PC = 2 h ) one, according to the estimations 
of the QM [H, [H^. Indeed, QM predictions show that 2 3 Pi and \ X D 2 states are the only ones compatible with 1 ++ 
or 2 \~ quantum numbers and lying approximately in the same energy region as the X(3872). The relativized QM 
predicts these states to be at energies of 3.95 and 3.84 GeV, respectively. However, the most recent results by the 
LHCb collaboration [lH seem to favor the 1 ++ quantum numbers. 

Our idea is thus to sec whether the introduction of loop corrections into the QM can help clarify the problem of 
the nature of the X(3872). Indeed, we think that the uncommon properties of the X(3872) are due to its proximity 
to the DD* decay threshold and cannot easily be explained within a standard quark-antiquark picture for mesons. 

In our calculation of Table [V] we have re- fitted the spectrum of charmonia through Eq. t)2b[) ; here, the mass of a 
meson results from the sum of a bare energy term computed within the relativized QM of Ref. [2[ , with a self energy 
correction computed within the unquenched quark model formalism of Refs. [3ll l50j. According to our results for the 
masses of the 2 3 Pi and l 1 !^ states, i.e. 3.908 and 3.741 GeV, respectively, the X(3872) is compatible with the meson 
Xci(2P) and includes an extra component due to the coupling to the meson- meson continuum, which is responsible 
for the downward energy shift. 

The second possibility is to treat the X(3872) as a DD* molecular state with 1 ++ quantum numbers [25l [29|. 
According to Refs. [73| . the DD* system with 1 ++ quantum numbers can be found by pion exchange and forms a 
meson molecule. More recent molecular model calculations [74| . including quark exchange kernels for the transitions 
DD* —> pJ/ty, wJ/\f r in order to predict the wJ/$ decay mode of the A(3872) [Hj], predict large isospin mixing due 
to the mass difference between D°D*° and D + D*~ . Nevertheless, in Ref. [25[ the authors observe that the one-pion 
exchange bin ding mechanism should be taken with greater caution in the DD* case than in the NN case (see also 
Refs. [Slzllzi)- 

Another important test for the properties of the X(3872) consists of estimating its strong and radiative decay rates 
[H, US, 53- In Ref. [4l|, the authors re-examine the re-scattering mechanism for the X(3872), which decays to J/i(ip(u>) 
through the exchange of D(*) mesons between intermediate states D(D) and D*(D*). Their results for the ratio 
R p /u> ~ 1= between the decay modes A(3872) -> J/ipp and A(3872) -> J/ipw, and for the rate A(3872) -> D°D°Tr°, 
favor a charmonium cc interpretation for the X(3872). In Ref. [io| . the author uses semi-quantitative methods 
to study some properties of the X(3872); he points out that the binding mechanism and the production rates are 
incompatible with the molecule interpretation. This is also suggested by the CDF II paper [77J where the authors 
observe also prompt production and discuss that a meson-meson molecule with a dimension of a few fm and intrinsic 
fragility cannot be prompt produced. By contrast, Refs. [25|, [2t| [3(| HI] suggest a molecular interpretation for the 
A(3872). 

Finally, we do not think that our arguments can, on their own, clarify the picture of the X(3872) resonance 
completely. Thus, it will be necessary to analyze other properties of this meson, such as strong and electromagnetic 
decays, in order to draw a definitive conclusion. In particular, we intend to calculate some of these observables within 
the UCQM [3l], by also taking the contribution of qq sea pairs into account [5l| . 

D. Discussion of the results 

In this paper we have presented the results of an unquenched quark model calculation of the self energy corrections 
to the spectrum of Iff, 25, IP, 2P and ID charmonium states. In the unquenched quark model, developed in the 
baryon sector in Refs. [53] and in the meson sector in Refs. [3l| , the effects of quark-antiquark sea pairs are introduced 
explicitly into the QM through a QCD-inspired 3 Po pair-creation mechanism. The UCQM model parameters are fitted 
to the reproduction of strong decay widths, as is shown in Sec. IIII Al 

The self energies are corrections to the bare meson masses arising from the coupling to the meson-meson continuum. 
Neglected in naive QM's, these loop effects provide an estimation of the quality of the quenched approximation used 
in QM calculations in which only valence quarks are taken into account. Something similar also happens in the case 
of lattice QCD, where one has to unquench the calculations in order to evaluate the contribution of the sea quarks to 
a certain observable. Therefore, one could say that these kinds of studies can be thought of as tests of the QM and of 
its range of applicability, and also as an enlargement of the model. Several studies on the goodness of the quenched 
approximation in the QM have already been conducted, such as those of Refs. [Ill, [Ho|, HH, in both the baryon and 
meson sectors. If the departure from the QM results is substantial, one can see new physics emerging or better extra 
degrees of freedom. This is the case of the X(3872), which in our picture can be described as a cc state plus higher 
Fock components mainly due to DD* and D*D* loops. 

Our results for the self energies of charmonia show that the pair-creation effects on the spectrum of heavy mesons 
are relatively small. Specifically for charmonium states, they are of the order of 2 — 6%, while we have shown in 
Refs. [3l| that the bottomonium mass shifts induced by the coupling to the meson-meson continuum are less than 
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approximately 1%. The relative mass shifts, i.e. the difference between the self energies of two meson states, are in 
the order of a few tens of MeV. However, as QM's can predict the meson masses with relatively high precision in the 
heavy quark sector - higher than can be obtained in the light meson sector or in baryon spectroscopy - even these 
corrections can become significant, such as in the case of the X(3872). 

It is interesting that the relative contribution of these corrections to meson masses decreases as the masses of the 
constituent quarks involved in the calculation increase. Moreover, qq pair creation is a relativistic effect, i.e. more 
important for low energy states. This is why we think that it would be quite interesting to use this formalism in the 
study of light mesons, for which relativistic effects, including qq pair creation, could make important corrections to 
the meson masses. 
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Appendix A: SUf(4) couplings 

The SUf(4) flavor couplings that we have to calculate in the 3 P model are (F B (14)Fc(32)|F A (12)F (34)) for the 
first diagram of Fig. [TJ and (Fb(32)Fc(14)|F<i(12)-Fo(34)} for the second diagram, where Fx{ij) represents the flavor 
wave function for the meson X (i.e. the initial meson A, the final mesons B and C or the 3 Pq created pair 0) made 
up of the quarks i and j. These overlaps can be easily calculated if we adopt a matrix representation of the mesons 
[20j. In this case, the two diagrams become, respectively, 

(F B (U)F C (32)\F A (12)F (M)} = Tr[F A FjF F£] 

= \Tr[F A F?FT] , 

(F B (32)Fc(14)|F A (12)F (34)) = Tr[F A F^F F^} [Al) 

= \Tr[F A F^Fl\ . 

For the SUf(5) flavor couplings, which have already been used for the bottomonium self energies in a preceding 
paper |3lj . the formulas arc 

(F B (14) J F 1 c (32)|F j4 (12)Fo(34)) - Tr [F A Fg FqFq] 

= -j=Tr{F A Fg F£] , 
(F B (32)F C (14)|F4(12)F (34)) = Tr[F A Fg F F?] 



(A2) 



±Tr[F A FZFT} 



Appendix B: Effective strength j 

It is known that the standard 3 Po model should not be applied for heavy-quark pair creation [3l| ; alternatively, 
the contribution from heavy channels should somehow be suppressed. Thus, in order to minimize the contributions 
from cc loops in Eq. (jfo]), we use the modified pair-creation mechanism of Refs. (23l. [3lj. This involves substituting 

the pair-creation strength of the 3 Pq model, 70, with an effective strength, 7^, defined as 

7o Cff =^7o, (Bl) 
rrn 

with i = n (i.e. word), s, c and b (see Table HT|) . 
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